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^ ■ Abstract 

A master equation for the deformed quantum harmonic oscillator interacting with 
a dissipative environment, in particular with a thermal bath, is derived in the micro- 
scopic model by using perturbation theory. The coefficients of the master equation 
and of equations of motion for observables depend on the deformation function. The 
steady state solution of the equation for the density matrix in the number represen- 
. tation is obtained and the equilibrium energy of the deformed harmonic oscillator is 

calculated in the approximation of small deformation. 
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S ■ 1 Introduction 

Q^. For more than a decade a constant interest has been induced to the study of deformations 
of Lie algebras - so-called quantum algebras or quantum groups, whose rich structure pro- 
^ ' duced important results and consequences in statistical mechanics, quantum field theory, 
5t , conformal field theory, quantum and nonlinear optics, nuclear and molecular physics. Their 

use in physics became intense with the introduction in 1989, by Biedenharn [1] and Mac- 
Farlane [2], of the g-deformed Heisenberg-Weyl algebra (g-deformed quantum harmonic 
oscillator). Since then the properties of the deformations of the harmonic oscillator have 
been investigated by many authors. Several kinds of generalized deformed oscillators have 
been introduced. There are, at least, two properties which make g-oscillators interesting 
objects for physics. The first is the fact that they naturally appear as the basic building 
blocks of completely integrable theories. The second concerns the connection between q- 
deformation and nonlinearity. In Refs. [3, 4, 5] it was shown that the g-oscillator leads to 
nonlinear vibrations with a special kind of the dependence of the frequency on the ampli- 
tude. For example, the g-deformed Bose distribution produces a correction to the Planck 
distribution formula [3, 4, 6]. 
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In the present paper we intend to study the connection of the quantum deformation and 
quantum dissipation, by setting a master equation for the deformed harmonic oscillator in 
the presence of a dissipative environment, which is shown to be the deformed version of 
the master equation obtained in the framework of the Lindblad theory for open quantum 
systems [7] . When the deformation becomes zero, we recover the Lindblad master equation 
for the damped harmonic oscillator [8, 9]. We are interested in describing the role of 
nonlinearities which appear in the master equation, this goal being motivated by the fact 
that the q'-oscillator can be considered as a physical system with a specific nonlinearity, 
called g-nonlinearity [3, 4]. For a certain choice of the environment coefficients, a master 
equation for the damped deformed oscillator has also been derived by Mancini [10]. In 
Ref. [11], EUinas used the g'-deformed oscillator for treatments of dissipation of a two-level 
atom and of a laser mode. 

The paper is organized as follows. In Sec. 2 we remind the basics about the generalized 
deformed quantum oscillator, in particular the /-oscillator and gr-oscillator. Using a variant 
of the Mancini's model [10], in Sec. 3 we derive a master equation for the /-deformed 
oscillator in the presence of a dissipative environment. The equations of motion obtained 
for different observables present a strong dependence on the deformation. Then in Sec. 
4 we write and solve in the stationary state the equation for the density matrix in the 
number representation. In the particular case when the environment is a thermal bath, we 
obtain an expression for the equilibrium energy of the oscillator in the approximation of a 
small deformation parameter. A summary and conclusions are given in Sec. 5. 

2 Deformed quantum oscillators 

It is known that the ordinary operators {1, a, a\ N} form the Lie algebra of the Heisenberg- 
Weyl group and the linear harmonic oscillator can be connected with the generators of the 
Heisenberg-Weyl Lie group. The generalized deformed quantum oscillators [12, 13] are 
defined as the algebra generated by the operators {1, A, A\ N} and the structure function 
F{N), which satisfy the following relations: 

[A,N]^A, [At,7V] = -At (1) 

and 

AA'' ^ F{N + 1), A^fA^FiN), (2) 

where F{N) is a positive analytic function with -F(O) = and N is the Hermitian number 
operator. It follows that the following commutation and anticommutation relations are 
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satisfied: 



[A,A^]=F{N + 1)-F{N), {A,A^} =F{N + 1) + F{N). (3) 

The structure function is a characteristics of the deformation. The number operator N is 
not equal to A^A as in the ordinary case. For F{N) — N one obtains the relations for the 
usual harmonic oscillator. Another choice is 

JV _ -N 

F{N) = _ = [N], (4) 

where the dimensionless c-number q is the deformation parameter. Then the operators 
A and are called g-deformed boson annihilation and creation operators [1, 2]. For the 
g-deformed harmonic oscillator the relations (2) become 

AA'^^[N + 1], A^A^[N], (5) 

with the commutation relation 

[A,A^] = [N+1]-[N]. (6) 

If q is real and positive, then 

^ sinh(iVlng) 
sinh(lng) 

and the condition of Hermitian conjugation [A^)"^ = A is satisfied. In addition to the 
commutation relation, there exists for the g-deformed oscillator the reordering relation 

AA^ - q'^^A^'A = q^^, (8) 

which is usually taken as the definition of g-oscillators. 

In the limit ? — > 1, g-operators tend to the ordinary operators because limg_+i[A'"] = N. 
Then Eqs. (6) and (8) go to the usual boson commutation relation [A, = 1. 

The g-deformed boson operators A and A^ can be expressed in terms of the usual boson 
operators a and (satisfying [a, a^] — 1^ N — a^a and [a, N] — a, [a^, N] — —a^) through 
the relations [14, 15]: 



V N+1 V V iV + 1 V 

Using a nonlinear map [16, 17], the g-oscillator has been interpreted [3, 4] as a non- 
linear oscillator with a special type of nonlinearity which classically corresponds to an 
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energy dependence of the oscillator frequency. Other nonlinearities can also be introduced 
by making the frequency to depend on other constants of motion, different from energy, 
through a deformation function / [3, 18]. Let us define the /-deformed oscillator operators 
[18] 

A = af{N) = f{N+l)a, = f{N)a^ = f{N + 1), (10) 

where N — a'a. They satisfy relations (1) and the commutation relation 

[A,A^] = {N + l)f{N+l) - Nf{N). (11) 

The function / has a dependence on the deformation parameter such that when the de- 
formation disappears, then / ^ 1 and the usual algebra is recovered. Without loss of 
generality, / can be chosen real and nonnegative and it is reasonable from the physical 
point of view to assume [5] that /(O) = 1 and f{N) = 1 for a suitable large N. A de- 
formation function depending on Laguerre polynomials has been used in [19]. A different 
type of deformation has been considered by Sudarshan [20] , who introduced the so-called 
harmonious states. Other examples of deformed oscillators are connected to the excited 
coherent states introduced by Agarwal [21] and Dodonov [22]. The transformation (10) 
from the operators a, a' to A, represents a nonlinear non-canonical transformation, since 
it does not preserve the commutation relation. The notion of /-oscillators generalizes the 
notion of ^-oscillators. Indeed, if 



sinh(iVlng) ^^^^ 



N \ A^sinh(lng)' 

then the operators A,A^ in Eqs. (10) satisfy the g-deformed commutation relations (6). 
This means that a Hamiltonian operator of the form A^A = f{N)a'^af{N) has a spectrum 
with the same structure as the spectrum of a^a. The difference is that the eigenvalues in the 
basis of the Fock space are np{n), n = 0,1, 2,..., instead of n. This spectrum associated 
with gf-deformation grows with n like sinh(nlng), i. e. exponentially for large occupation 
numbers n, compared to the ordinary case, in which the spectrum is equidistant. The 
Hamiltonian of the /-deformed harmonic oscillator is {uj is the ordinary frequency) 

n = ^{AA^ + A^A) = ^[{N + l)f{N + 1) + Nf\N)\. (13) 
It is diagonal on the eigenstates \n > and in the Fock space its eigenvalues are 

En = ^[{n + l)f{n + 1) + nf{n)]. (14) 
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In the limit / — > 1 (g — > 1 for g-oscillators) , we recover the ordinary expression En — 
nw{n + 1/2). 

Using the operator Heisenberg equation with the Hamiltonian (13) 

in^^-[a,n] (15) 

or the evolution operator U{t) — exp[—{i/h)H{N)t\, we obtain the following solutions to 
the Heisenberg equations of motion for the operators a and [10, 18]: 

a(t) — exp[—iu!fl{N)t]a, a' {t) ^ exp[iujn{N)t], (16) 

where 

n{N) ^1[{N + 2)f{N + 2) - Nf{N)]. (17) 
For a g-deformed harmonic oscillator, 

m)-li[N + 2]-[N]) (18) 
and for a small deformation parameter r (r = In g) , 

n{N) = l + ^{N + lf. (19) 

3 Quantum Markovian master equation 

In order to discuss the dynamics of the open system S, we use a microscopic description 
of the composite system S+B. As the subsystem S of interest we take the /-deformed 
harmonic oscillator with the Hamiltonian (13), and B is the environment (bath) with 
the Hamiltonian Hb- The coupled system with the total Hamiltonian Ht — 'H-\- Hb + V 
{V is the interaction Hamiltonian) is described by a density operator x(i), which evolves 
in time according to the von Neumann-Liouville equation 

^^4|H..XW1. (20) 

When the Hamiltonian evolution of the total system is projected onto the space of the 
harmonic oscillator, the reduced density operator of the subsystem is given by 

p{t)^Tr^x{t). (21) 

The derivation of the reduced density operator in which the operators of the environment 
system have been eliminated up to second order of the perturbation theory can be taken 
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from literature [23, 24, 25, 26, 27]. We assume that the interaction potential V is linear in 
the coordinate operator si — q and momentum operator S2 — p in the Hilbert space of the 
subsystem. Then, following [26], we can write down the master equation for the density 
operator of the open quantum system in the Born-Markov approximation: 

^ = ~[n,p{t)] + ^ E //t'{c*.(Oh,pW^.(-0] + ^^^.(OK(-^')pW,^.]},(22) 

where the coefficients Cij{t) are correlation functions of the environment operators. It 
is assumed that the correlation functions decay very rapidly on the time scale on which 
p{t) varies. Ideally, we might take Cij{t') ~ S{t'). The Markov approximation relies on 
the existence of two widely separated time scales: a slow time scale for the dynamics of 
the system S and a fast time scale characterizing the decay of environment correlation 
functions [26]. 

In order to get the time dependence of the operators Si{t) = q{t) and S2{t) = p{t), we 
express them through the relations (m is the oscillator mass) 



= V2^^"^^^^ + "^^^^' P^^^ = '\l^^''^(^^-<^^^ (23) 
and then insert Eq. (16) for a{t) and a^(i). Then the master equation results 

+ /* dt'{CUt')[q, p{t){qE^ + E+q - —{pE_ - E^p))] 

+Ci,{t')[{qE_ + E+q - :^{pE- - E+p))p{t), q] 
+iC;^{t')\p, p{t){mu{qE^ - E+q) - i{pE_ + E+p))] 
+iC22{t')[{mcu{qE_ - E+q) - i{pE^ + E+p))p{t),p] 
+iC*^it')[q, p{t){ww{qE_ - E+q) - i{pE_ + E^p))] 

+C2i{t')[{qE_ + E+q - —{pE_ - E+p))p{t),p\ 



+CUt')\p, p{t){qE_ + E+q - —{pE- - E+p))] 
+iC,2it')[imwiqE_ - E+q) - i{pE_ + E+p))p{t), q]}, (24) 

where we have introduced the following notations: 

E+ = exp[iujQ{N)t'], E_ = exp[-iu;Q{N)t']. (25) 

If the environment is sufficiently large we may assume that the time correlation func- 
tions decay fast enough to zero for times t' longer than the relaxation time of the 
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environment: t' ^ ts- Therefore, if we are interested in the dynamics of the subsystem 
over times which are longer than the environment relaxation time, t ^ ^b, we may use the 
Markov approximation and replace the upper limit of integration f by oo. Physically, this 
amounts to assuming that the memory functions Cij{t')E±{t') decay over a time which is 
much shorter than the characteristic evolution time of the system of interest. After certain 
assumptions [23, 27], one can define the complex decay rates, which govern the rate of 
relaxation of the system density operator as follows: 

dt'Cu{t')E+ = dt'C*,,{t')E+ = Dpp{Q), (26) 

fOO POO 

dt'C22it')E+ = I dt'C*,it')E+ = D,,iQ), (27) 

/ dt'Cu(t')E+ = / dt'C;,{t')E+ = -Dp,(n) + -A(Q), (28) 
Jo Jo Z 

with Dpp{n) > 0, Dqq{n) > and 

D,,{n)D,,{n) - Dl^{n) > ^x'{n). (29) 

In fact, Dpp{Q), Dqq{Vi), Dpq{Vt) and play the role of deformed diffusion and, respec- 
tively, dissipation coefficients and the relation (29) ensures the positivity of the density 
operator. The existence of these coefficients reflects the fact that, due to the interaction, 
the energy of the system is dissipated into the environment, but noise arises also (in par- 
ticular, thermal noise), since the environment also distributes some of its energy back to 
the system. In addition, we assume in the following X{Q) — \ — const. Then the master 
equation (24) for the damped deformed harmonic oscillator takes the form 

dp 1 

+ ^{[{{Dpp{^). q] + —[Dpp{n),p\)p, q] + [({Dg,(0),p} - imw[Dqq{Q), q])p,p\ 

+[{rmj[iDpq{rt) + ^A, q] - {Dpg{n) - yA,p})p, q] 

-[{ — [zDpq{n) - h,p] + {Dpq{n) + ^A, q})p,p] + H.c.}. (30) 

We notice that the deformation is present in both the commutator containing the oscillator 
Hamiltonian 7i, as well as in the dissipative part of the master equation, which describes 
the influence of the environment on the deformed oscillator. This master equation preserves 
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the Hermiticity property of the density operator and the normahzation (unit trace) at all 
times, if at the initial time it has these properties. In the limit / — > 1 (Q — > 1), the 
deformation disappears and Eq. (30) becomes the Markovian master equation for the 
damped harmonic oscillator, obtained in the Lindblad theory for open quantum systems, 
based on completely positive dynamical semigroups [7, 8, 9]. 

Expressing the coordinate and momentum operators back in terms of the creation and 
annihilation operators and introducing the notations 

D^{n)^^[mujD,,{n) + ^^], D_{n) ^ hmujD,,in) - ^^], (31) 
the master equation (30) for the damped deformed harmonic oscillator takes the form: 

+{[[D+(0)a,p],at] - [[a\D_{n) + ^^,,(0)), p], a^] - ^[a\{a, p}] + H.c.}. (32) 

Mancini considered in Ref. [10] a squeezed bath for the dynamics of the damped deformed 
harmonic oscillator and his model can be recovered by taking the following coefficients in 
Eq. (32): 

L'+(Q)=7(A^+^), D_{Q)+'-D,,{Q)^-^M, A = 7- (33) 

In the particular case of a thermal equilibrium of the bath at temperature T {k is the 
Boltzmann constant), we take the diffusion coefficients of the form (in concordance with 
Mancini's results [10]) 

mwD,,{n) = = h coth D,,{n) = 0. (34) 

In the limit Q — > 1, the deformed diffusion coefficients (34) take the known form obtained 
for the damped harmonic oscillator in the particular case when the asymptotic state is a 
Gibbs state [8, 9]: 

hmuj ^ hoj ^ Ti . ^ Tiuj ^ ^ 
J^vv — A coth — — , Daa — A coth — — , D^a — 0. (35) 

The meaning of the master equation becomes clear when we transform it into equations 
satisfied by the expectation values of observables involved in the master equation, < O >— 
Tr[p[t)0], where O is the operator corresponding to such an observable. Wc give an 
example, multiplying both sides of Eq. (32) by the number operator and taking the 
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trace. In the case of a thermal bath, with the diffusion coefficients (34), the equation of 
motion for the expectation value of N has the form 

I < iV >= A[< (coth ^ _ l)(iV + 1) > - < (coth + i)jv ». (36) 

This equation leads to a time dependence of the number of quanta on dissipation and 
temperature, compared to the case of an oscillator without dissipation, where the number 
of quanta is conserved. We remark that in the case of a thermal bath at T = 0, Eq. (36) 
takes the form 

j^<N >^-2X<N >, (37) 

so that the average number of quanta < N > does not depend on deformation, it only 
decreases exponentially with dissipation. 

We consider another example, taking the simplest case of a thermal bath at T = 0, 
when both diffusion and dissipation coefficients do not depend on the deformation, = 
A/2 = const. Even in this situation, the equations of motion for the expectation values are 
yet complicated, because they do not form a closed system. Multiplying both sides of Eq. 
(32) by the operators a and, respectively, Q{N)a and taking throughout the trace, we get 
the following equations for the expectation values of these operators: 

< a >= -luj < Q{N)a > -\ < a >, (38) 

< n(N)a >= -ioo < n^(N)a > +X < [2Nn(N - 1) - (2A^ + l)n(N)]a > . (39) 

These examples show that the equations of motion contain nonlinearities introduced by 
the deformed Hamiltonian Ti and, therefore, depend on the deformation function. 

4 Equation for the density matrix of the damped de- 
formed oscillator 

Let us rewrite the master equation (32) for the density matrix by means of the number 
representation. Specifically, we take the matrix elements of each term between different 
number states denoted by \n >, and using N\n >— n\n >, a'^\n >— \Jn + \\n + 1 > and 
a\n >— \/n\n — 1 >, we get 
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^ = -y [m/2(m) + (m + l)fim + 1) - nf{n) - {n + l)fin + 
-[(m + l)D+{n{m)) + mL'+(fi(m - 1)) + (n + l)L)+(r](n)) + nD+(fi(n - 1)) - A]p^„ 

+^(m + l)(n+l)[L'+(0(m)) + D+{n{n)) + X]pm+i,n+i 
+V^[D+{n{m - 1)) + D+{n{n - 1)) - X]pm-i,n-i 
-^J{m+l)n[D_{n{m)) + D_((](n - 1)) - ^{Dpg{nim)) + Dpg{n{n - l)))]pm+i,n-i 

m— l,n+l 

+ ^{m + l){m + 2)[D_{Q{m + 1)) - ^/^p,(i^(m + l))]pm+2,n 
+yJ{n + l){n + 2)[D4n{n + 1)) + ^Dpg{n{n + l))]p 

m,n+2 

+ ^mim-l)[D_{n{m - 2)) + ^Dp,{n{m - 2))]p^_2,„ 

+^n{n-l)[D4n{n - 2)) - ^Dp,{n{n - 2))]p„,„_2.(40) 

Here, we have used the abbreviated notation pmn =< iTi\p{t)\n > . This equation gives an 
infinite hierarchy of coupled equations for the matrix elements. When 

L'_(fi(n)) = 0, Dpg{Q{n))^0, (41) 

the diagonal elements are coupled only amongst themselves and not coupled to the off- 
diagonal elements. In this case the diagonal elements (populations) satisfy a simpler set 
of master equations: 

^^M. = _[2(n + l)D+{Q{n)) + 2nD+{n{n - 1)) - \]P{n) 
+{n + l)[2D+{n{n)) + X]P{n + 1) + n[2D+{n{n - 1)) - X]P{n - 1), (42) 

where we have set P{n) = Pnn- We define the transition probabilities 

t+{n) = {n + l)[2D+{n{n)) - X], t-{n) = n[2D+{n{n - 1)) + X]. (43) 

With these notations Eq. (42) becomes: 

= t+(n - l)P(n - 1) + t4n + l)P(n + 1) - [t+{n) + i_(n)]P(n). (44) 
The steady state solution of Eq. (44) is found to be 
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We note that in the steady state the detailed balance condition holds: 



t-{n)P{n) =t+{n-l)P{n-l). (46) 

In the particular case of a thermal state, when the diffusion coefficients have the form 
(34), the stationary solution of Eq. (44) takes the following form: 

= Vexp{-^[(n + l)f (n + l)+nf (n)]}, (47) 

where 

Zr = P(0)e./^ (48) 
and Zf is the partition function: 

oo t . 

Zf = E exp{-^[(n + l)f{n + 1) + nf{n)]}. (49) 

In his model, Mancini obtained a result similar with Eq. (47). 
Using Eq. (14), the distribution (47) can be written 

P'ssin) = Zj'exp{-^). (50) 

Expression (47) represents the Boltzmann distribution for the deformed harmonic oscilla- 
tor. In the limit / — > 1 the probability Pgg{n) becomes the Boltzmann distribution for the 
ordinary harmonic oscillator with the well-known partition function 

(51) 



2sinh^' 



For the g-oscillator described by the Hamiltonian H (13) and weakly coupled to a 
reservoir kept at the temperature T, the gr-deformed partition function can be obtained 
as a particular case of the partition function Zf (49), by taking the deformation function 
(12): 

^ huj sinh(r(n + 1)) + sinh(rn) 

In the limit of a small deformation r we can write [4] Zg = Z + br^ , where 

6 = -^(2n3 + 3n2 + n), (53) 
1 2 
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with 

1 e^ + 1 e2/^ + 4e^ + l ^ hu 

^=^' "^F^' " = (e/^-l)3 ' ^ = ^^^^ 

We can calculate the equilibrium energy by using the formula 

1 f)7 



and obtain 



where 



^(oo) = — (coth^ + r^c), (56) 



- (^^T^^vT-]- (57) 



(e/3-l)2Le/3_i ^ (6/5-1)2 

We note that in the approximation of a small deformation parameter r, the energy of the 
deformed damped oscillator depends on the oscillator ground state energy hu}/2 and on the 
temperature T. Evidently, when there is no deformation (r 0), one recovers the energy 
of the ordinary harmonic oscillator in a thermal bath [8, 9]. In the limit T" — > 0, one has 
c — > 0, E{oo) — huj/2 and the deformation does not play any role. 

5 Summary and conclusions 

Our purpose was to study the dynamics of the deformed quantum harmonic oscillator in 
interaction with a dissipative environment, in particular with a thermal bath. We derived 
in the Born-Markov approximation a master equation for the reduced density operator of 
the damped /-deformed oscillator. The one-dimensional /- or g-oscillator is a nonlinear 
quantum oscillator with a specific nonlinearity and, consequently, the diffusion and dissi- 
pation coefficients which model the influence of the environment on the deformed oscillator 
depend strongly on the introduced nonlinearities. The equations of motion for the observ- 
ables of the considered system are also nonlinear. In the limit of zero deformation, the 
master equation for the deformed damped oscillator takes the form of a master equation for 
the damped oscillator obtained in the framework of the Lindblad theory of open quantum 
systems based on quantum dynamical semigroups. We have also derived the equation for 
the density matrix in the number representation. In the case of a thermal bath we obtained 
the stationary solution, which is the Boltzmann distribution for the deformed harmonic 
oscillator. In the approximation of a small deformation, we obtained the expression of the 
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equilibrium energy of the deformed harmonic oscillator, which depends on the oscillator 
ground state energy and on the temperature. 

The master equation for the damped deformed harmonic oscillator is an operator equa- 
tion. It could be useful to study its consequences for the density operator by transforming 
this equation into more familiar forms, such as partial differential equations of Fokker- 
Planck type for the Glauber, antinormal ordering and Wigner quasiprobability distribu- 
tions or for analogous deformed quasiprobabilities [5] associated with the density operator. 
It could also be interesting to find the states which minimize the rate of entropy produc- 
tion for the damped deformed harmonic oscillator. In the case of the undeformed damped 
oscillator such states are represented by correlated coherent states. For the damped de- 
formed oscillator the corresponding states could be deformed (nonhnear) coherent states 
[18] , which may play an important role in the description of the phenomenon of environment 
induced decoherence. The dissipative dynamics of deformed coherent states superposition 
and the related coherence properties have been studied recently by Mancini and Man'ko 
[28]. 
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